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Abstract 

In this paper, we employ the bifurcation theory of planar dynamical systems 
to investigate the travelling-wave solutions to a dual equation of the Kaup- 
Boussinesq system. The expressions for smooth solitary-wave solutions are 
obtained. 
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1. Introduction 

Since the theory of solitons has very wide applications in fluid dynamics, 
nonlinear optics, biochemistry, microbiology, physics and many other fields, 
the study of soliton solutions to nonlinear partial differential equations has 



It is known that, 



become an increasingly important area of research 
solitons are the solitary waves that retain their individuality under interaction 
and eventually travel with their original shapes and speeds. Therefore, to 
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investigate the soliton solutions, one must firstly obtain the solitary-wave 
solutions. Many efforts have been denoted to seekingsolitary-wave solutions 
to nonlinear partial differential equations (see, e.g., Q]- lsj]). 

Recently, Guha 1 14] studied the dual counter part of the following Kaup- 
Boussinesq system |l5fl . 



u t = v xxx + 2(uv) x , 
v t = u x + 2vv x , 



(1-1) 



where u(x, t) denotes the height of the water surface above a horizontal bot- 
tom and v(x,t) is related to the horizontal velocity field. Using moment of 
inertia operators method and the frozen Lie-Poisson structure, Guha derived 
the dual counter part of system (jl.ip . that is 



m t + (mu + \u\ + u 2 + 2uv) 3 
Pt + (pu)x = 0, 



0, 



;i.2) 



where m 



u x . 



+ fxu + Xv, p = Xu + v. System (11.21) is a two component 



integrable system [141 ] . When fi = X = 0, it becomes 



Uxxt + {uu xx + \u\ + u 2 + 2uv) 3 
v t + (uv) x = 0. 



;i.3) 



Various aspects of the Kaup-Boussinesq system (11.11) have been studied. 



For instance, Smirnov 



lq| obtained real finite gap regular solutions to system 



(11. II) . and Borisov et al. [17j applied the proliferation scheme to system ( II. ip . 
Also the closely related variant of system (11.11) have been studied intensively 



see 



18|-[29|]). However, it seems that the dual equation of system (11.11) has 



attracted little attention. 
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In t 

(see 



lis paper, we use the bifurcation theory of planar dynamical systems 



30|-|33J]) to investigate the travelling-wave solutions to system (jl.3p and 
obtain analytic expressions for its smooth solitary-wave solutions. To the 
best of our knowledge, the solitary-wave solutions to system (11.31) have not 
been reported in the literature. The bifurcation method was first used by Li 
and Liu [34J to obtain smooth and non-smooth travelling-wave solutions to a 
nonlinearly dispersive equation and was later employed by many authors to 
derive a variety of travelling-wave solutions to a large number of nonlinear 



partial differential equations 



35|-|42| 



The remainder of the paper is organized as follows. In Section 2, using the 
travelling- wave ansatz, we transform system (11.31) into a planar dynamical 
system and then discuss bifurcations of phase portraits of this system. In 
Section 3, we obtain the expressions for smooth solitary-wave solutions to 
system (11.31) . A short conclusion is given in Section 4. 

2. Bifurcation analysis 

Let £ = x — ct, where c is the wave speed. By using the travelling-wave 
ansatz u(x, t) = <p(x — ct) = <p(£), v(x,t) = ip(x — ct) = we reduce 

system (11.31) to the following ordinary differential equations: 



-off + w + y 2 + ^ + i<p*i>y = 0, 

-ap' + (tpipy = 0. 

Integrating (12.11) once with respect to £, we have 

-off' + tpip" + y 2 + tp 2 + 2tpip = gi , 
-cip + <pip = g 2 , 



(2.1) 



(2.2) 



where g\,Q2 are two integral constants. 

From the second equation in system (I2.2p . we obtain 

^ = ^^. (2.3) 
ip — c 



Substituting (12. 3p into the first equation in system (12. 2p yields 

^ = {p - c){ 9l - ±<p> 2 - <p 2 ) - 2g 2 ^ (2 4) 



Let p' = ^y, then we get the following planar dynamical system: 
dp - V2.. 

dy _ V2{{p - c){ 9l - \y 2 - p 2 ) - 2g 2 p) ( 2 - 5 ) 

This is a planar Hamiltonian system with Hamiltonian function 

H{tp, y) = p 4 - j p 3 + (4(7 2 - 2 9l )p 2 + Acgup + (p - cfy 2 = h, (2.6) 

where h is a constant. 

Note that (12.51) has a singular line (p — c. To avoid the line temporarily we 
make transformation d£ = — c) 2 d(. Under this transformation, Eq. (12.51) 
becomes 

a > (2.7) 



2' 

System (12.51) and system (12.71) have the same first integral as (12.61) . Con- 
sequently, system (12.71) has the same topological phase portraits as system 
(12.51) except for the straight line tp — c. 

For a fixed h, (12.61) determines a set of invariant curves of system (12. 7p . As 
h is varied, (I2.6p determines different families of orbits of system (12 .7p having 



different dynamical behaviors. Let M((p e , y e ) be the coefficient matrix of the 
linearized version of system (\2.7\i at the equilibrium point ((p e ,y e ), then 



(ip e -c)y e Uy e -c) 
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M{<p e ,y e )= | vr - ' 2 | (2.8) 

-3ipl + 2ap e + gi - 2g 2 - ^y 2 -(ip e - c)y e 



and at this equilibrium point, we have 

J(<f e ,ye) = det M((p e ,y e ) 



■{(p e - cfy 2 e + ~(<p e - c) 2 (3^ - 2cp e -gi + 2g 2 + ^y 2 ) 



(2.9) 



p((p e , y e ) = trace(M(v? e , y e )) = 0. (2.10) 

It is easy to see that the equilibrium point of system (12.71) is in the form 
of (tp e ,0). At this equilibrium point, we have J((p e ,0) = \{ip e — c) 2 (3ipl — 
2cif e — g± + 2g 2 ) . By using the bifurcation theory of planar dynamical system, 
we know that if J(<p e , 0) > (or < 0), then the equilibrium ((p e , 0) is a center 
(or saddle) point; if J(cp e ,0) = 0, and the Poincare index of the equilibrium 
point is 0, then it is a cusp. 

Usually, a solitary-wave solution to system (11. 3p corresponds to a homo- 
clinic orbit of system (12. 71) . Therefore, to obtain solitary- wave solutions to 
system (11.31) . we need only to seek homoclinic orbits of system (12.71) and so 
only the saddle points are of interest. Firstly, we need to look for the possible 
zeros of the function 

f(<p) = -ip 3 + ap 2 + cup + b, (2.11) 

where a = g\ — 2g 2 , b = —cg\. 

In order to find all possible zeros of f{<p), we set 

f'(ip) = Sip 2 + 2ap + a = 0. (2.12) 



When A = 4c 2 + 12a > 0, we find two real roots to Eq. fl2.12j) as follows: 



<Pi = o( c ~ Vc 2 + 3a), (2.13) 



<p* 2 = -(c + Vd* + 3a), (2.14) 

with <p\ < p* 2 . When A = 4c 2 + 12a < 0, the inequality f(<p) < holds. In 
this case, if (<p e , 0) is an equilibrium point of system (12. 71) . then it is a center 
point (or a cusp) because J(tp e ,0) = — \(<p e ~ c ) 2 f{<Pe) > 0. Therefore, in 
the following, we always suppose A = 4c 2 + 12a > 0. 

Substitute fl2TT31 and f!2TT4j) into (|2T[jl . respectively, we get 

h = m) = %- 2 ^v^tm- 2 -^v^+m+j + b, (2.i5) 



h = m) = + ^ v^+3^ + f v^+3^ + y + 6, (2.16) 



with 



h ~ h = -|^(c 2 + 3a) < 0. (2.17) 
The equilibrium points of system (12. 7p have the following properties. 
Theorem 2.1. (%) // / 2 < 0, then system ( 2.7 ) has only one equilibrium 



point, denoted by (ipi,0)(ipi < ip* < tpl^), which is a center point; 



(2) If fi > 0, then system (2.1 ) has only one equilibrium point, denoted 
by (<f2, 0)(<£>* < p\ < which is a center point; 

(3) If fi = 0, then system \2. J\j has two equilibrium points, denoted by 
(v?3,0), (v?4,0)(c/?3 = tp* < ip* < y> 4 ). (ps,0) is a cusp, while (^4,0) is a center 
point; 
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(4) If fi = 0, then system \2. 7| ) /ios two equilibrium points, denoted by 
(<£>5,0), (y9 6 ,0)(^5 < <y?* < Lp\ = <f 6 ). ((p 5 , 0) a center point, while (y?6,0) zs 
a citsp; 

fjj If f\ < 0, / 2 > 0, i/ien system {2. 7| ) /ias t/iree equilibrium points, 
denoted by (<p 7 , 0), (<^ 8 , 0) and (<^9,0)(<^ 7 < 97* < y? 8 < ^ < y?g). (</?7, 0) and 
(</?9, 0) are two center points, while (ip$, 0) is a saddle point. 

In this paper, we only consider the case c > because in the case c < 
we will get analogous result. In order to give the details of the bifurcation, 
we fix the parameter a = — 1. Thus, we obtain the following two bifurcation 
curves of system (12.71) . 

2c 3 2c 2 



The above bifurcation curves divide the parameter space into three re- 
gions (see Fig. [T]) in which different phase portraits exist. By theorem 12.14 
we can see that only in regions (II), can system (12. 7p has saddle points. See 
Fig. [2] for an example of the corresponding phase portraits. 

3. Solitary-wave solutions to system ( 11.31) 

From the discussions in Section 2, we can see that, when the parameters 
a = —1, (6, c) G (II), system (12. 7\i has infinite many saddle points. So 
there are infinite many homoclinic orbits and system (II. 3p has infinite many 
solitary-wave solutions accordingly. 
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Figure 1: The bifurcation sets and bifurcation curves of system ( 12.7ft for the 
parameter c > \/3. 



In order to obtain exact expressions for solitary-wave solutions, we fix 

27 ^ 3 

Case I: V3 < c < 3 

In this case, there are two homoclinic orbits connecting with the saddle 
point (|,0), see Fig. 0(a) for an example. The two homoclinic orbits of 
system (12.71) or (12.51) can be expressed respectively as 



(v-f)v-^ + ^+¥-2 



if — c 



for (p_ < ip < -, 



(3.1) 



-^ + f^+f-2 

for - < (p < ip + , 



ip — c 



(3.2) 



where cp± = |(c ± V6c 2 — 18). 

Substituting (I3.ip . (13.21) into the first equation of system (12.51) . respec- 
tively, and integrating along the corresponding homoclinic orbit, we have 



s — c 



'v- (s - 



-s 2 + ^s+ 5 4-2 



(3.3) 
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(a) (b) 

Figure 2: The phase portraits of system (12.71) when the parameters c > \/3, 
a = -1 and b = -f^ + f . (a) c = 1.8; (b) c = 4. 



^ s-c 



^ = ~\H (3-4) 



«> (s-f)i/-s 2 + fs + ^-2 2 
It follows from (ET3)) . f[3T4j) that 

^ + arctan(a(y))+ ^ ln(/%,)) = ffifl, y_ < ^ < |, (3.5) 



and 



7r 2c \/2 c 

--arctan(a(^))--^====ln(-/?(^)) = — - < <p < <p+, (3.6) 



where 



a(<p) = ^ C (3.7) 

' -9p 2 + 6cy? + 5c 2 - 18 



_, , V6C 2 - 18 + ^-9^ + 6c^ + 5c 2 - 18 

/%) = ^ • (3-8) 

c — 099 

Therefore, we obtain two solitary-wave solutions to system ( 11.31) in the 
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following parametric forms: 

£ = ±V2(f + arctan(a(^)) + 7 ^=ln(^)), / ^ ^ c 

(p = <p, 



(V-<V< 3), (3-9) 



e = ±v / 2(| + arctan(a(^)) + ^=ln(/3(^)), c 

2 V6c ~ 18 < y» < -), (3.10) 

ib = - 33 -, 3 

r </3 — C ' 

and 

^i^f-arctan^-^^M-^)), c 

v (- < ^ < (3.11) 

p = <p, 

£ = ±V2($ - arctan(a(y?)) - 2c = ]n(-8(ip)), c 

= _92_ 3 

Now we take a set of data and employ Maple to display the graphs of the 
above obtained solitary-wave solutions in Fig. [3J 
Case II: c > 3 

In this case, there is one homoclinic orbit connecting with the saddle 
point (§,0), see Fig. H^b) for an example. 

Similar to the Case I, we can obtain a solitary-wave solution to system 
(11.31) . given as (13. 9p . (13.101) . A typical such solution is shown in Fig. HI 

4. Conclusion 

In summary, by using the bifurcation method, we obtain analytic ex- 
pressions for smooth solitary-wave solutions to a dual equation of the Kaup- 
Boussinesq system (11.31) . The results of this paper suggest that, in addition to 
solving many single-component partial differential equations, the bifurcation 
method can be used to obtain travelling-wave solutions of two-component 
systems. 

10 





0.6- 
\ 0.5- 






\o.4- 






0y3- 





-30 -20 -10 10 20 30 

(a) 




(c) 



-30 -20 -10 10 20 30 



-0.66H 




-oho-_ 




¥ / : 




70.75- 




y-o.8o^ 





(b) 



-20 -10 



(d) 



10 20 



I.I. 

^^X-0.9- 


1 1 


\.o- 




V \ - 




-1U- 




-1.3r 





Figure 3: Solitary- wave solutions to system (11.31) when the parameters c 
1.8, a = -1, b = 0.168000 and g 2 = 1. 
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